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Several convenient formulae for the entanglement of two indistinguishable delocalised spin-i par- 
ticles are introduced. This generalizes the standard formula for concurrence, valid only in the limit 
of localised or distinguishable particles. Several illustrative examples are given. 
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Entanglement is a well-defined quantity for two dis- 
tinguishable qubits in a nonfactorizable quantum state, 
where it may be uniquely defined through von Neuman 
entropy and concurrence QIMIMQ- However, amongst 
the realistic systems of major physical interest, electron- 
qubits have the potential for a much richer variety of en- 
tanglement measure choices due to both their charge and 
spin degrees of freedom. For example, in lattice fermion 
models such as the Hubbard dimer, entanglement is sen- 
sitive to the interplay between charge hopping and the 
avoidance of double occupancy due to Hubbard repul- 
sion, which results in an effective Heisenberg interaction 
between adjacent spins In systems of identical par- 
ticles the main challenge is to define an appropriate en- 
tanglement measure which adequately deals with multi- 
ple occupancy states 0, S S ^> 0| ■ In the case of 
fcrmions such a measure must also account for the effect 
of exchange |0 as well as of mutual electron repulsion. 

Entangled fermionic qubits can be created with 
electron-hole pairs in a Fermi sea |0| and in the scat- 
tering of two distinguishable particles 0|. A spin- 
independent scheme for detecting orbital entanglement 
of two-quasiparticle excitations of a mesoscopic normal- 
superconductor system was also proposed recently . 

A consensus regarding the appropriate generalization 
of entanglement measure which would consider spin and 
orbital entanglement of electrons on the same footing 
has not, however, been reached yet. In any realistic 
solid-state device, spin entanglement is intimately related 
to the orbital degrees of freedom of the carriers, which 
cannot be ignored, even in otherwise pure spin entan- 
glement observations. In this paper we introduce spin- 
entanglement measure formulae valid for real electrons 
and show how, in general, spin-entanglement depends in 
an essential way on spatially delocalised orbitals. 

For two distinguishable particles A and B, each de- 
scribed with single spin-i (or pseudo spin) states s 
or I and in a pure state = J2ss' "ss' I^UIsO-B con- 

currence as a measure of entanglement is given by 



C = 2|a||aj4. - a|j.axi|. (1) 
Concurrence is related to the density matrix of a pair 



of spins 1^ and can be expressed in terms of spin- 
spin correlators (4'ab|'S'^<S'^|^'ab) and expectation val- 
ues (*ab|S'^(b)I*ab), where S^^^^ for A = x,y,z are 
spin operators corresponding to spin ^ or i?, respectively. 
This approach has proved to be efficient in the analysis 
of entanglement in various spin-chain systems with inter- 
action HIEIEIEI- 

Consider now the general problem of two interact- 
ing electrons in a pure state. It is clear that in some 
circumstances this system reduces approximately to an 
equivalent system of two interacting spins, for which the 
above entanglement formula is appropriate. Further- 
more, in the general case, entanglement between the spins 
of the fermions relates to measurements of spin irrespec- 
tive of their orbital motion. We consider therefore spin- 
entanglement for a general class of two-electron states on 
a lattice of the form 
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where c|^^ creates an electron with spin s on site i and N is 
the total number of sites. The system in question could 
be, for example, a tight-binding lattice containing two 
valence electrons occupying non-degenerate atomic or- 
bitals, or two electrons in the conduction band of a semi- 
conductor, for which the sites represent finite-difference 
grid points. In cither case, the interaction between the 
electrons is included together with any externally applied 
potential. 

The two electrons are in separate regions of space 
(measurement domains) [A] and [B] as illustrated in 
Fig.^a). Entanglement might be produced, for example, 
when two initially unentangled electrons in wave pack- 
ets approach each other and interact [Fig.^b)] and then 
again become well separated into distinct regions [A] and 
[B] [Fig. ^c)]. Here one should realise that in real mea- 
surements of entanglement, indistinguishable electrons 
would be detected and the formalism relevant to dis- 
tinguishable spins is not directly applicable. Neverthe- 
less, complete information regarding the spin properties 
of such a fermionic system is contained in spin correlation 
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functions for the two domains. The spin-measuring ap- 
paratus would measure spin correlation functions for two 
domains [A] and [B] rather than for two distinguishable 
spins A and B. 

Concurrence as a measure of entanglement for two elec- 
trons is related to the eigenvalues of the non-Hermitian 
matrix pp, where p is reduced density matrix given 
in terms of the electron spin correlations correspond- 
ing to the domains, and p is the time- reversed den- 
sity matrix as in Ref. 0. In general the eigenvalues of 
pp can be determined only numerically and a closed 
form for concurrence can not be obtained, unless the 
system exhibits additional symmetries. Possible sym- 
metries are conveniently studied through spin-spin cor- 
relation functions. We express spin operators for do- 
mains [A] and [B] with fermionic operators as the sum 
of operators for sites i within the domain [A] (or [B]), 
i.e., S\ = 5E^e[A]Ess'4^s.'C^.', where are Pauh 
matrices. For axially symmetric problems 20], where 



i^Kjnfl^) = and (vl^l^^.^^^.^^^^ 
rence may be written as 



= 0, concur- 



max(0,C|x,q|), 

2\{S+S^)\-2J{Plpl){Pipj,) 



C 

q, = 2\{S+S+)\-2J{P\P^){PXPl,), 



where S 
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are spm 



raising operators for domains [A\ or [B] and P 
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Y.i£A(B) ~ with n.is = c\^c^s, are spin-s pro- 

jectors operating in domains [A\ (or [B]). Fermionic ex- 
pectation values required in Eq. ij^J are then given in 
terms of the amplitudes in the normalised 1 5') as 
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(4) 



where the summation in Eq. |0J extends over all pairs 
[ij] such that i G [A] and j e [B]. In analogy to the 
Bell basis Q one can introduce (pf, = (V'jj ± ip]\)/V2 



and xfj = Wj ± V'ij)/V2 . ipfj, e.g., are the ampli- 



tudes for creating two electrons in a delocalised singlet 
or triplet state with zero total spin projection. It then 
follows from Eqs. l(21and01) that the electrons are com- 
pletely entangled, when either(i) ip^j — lapfj, x% = 0, or 

(ii) xTj — ^'^xtj^ ^ij ~ where c is a real constant. In 
the general case (i.e., without spin symmetries) C = 1 if 
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Figure 1: (Color online) (a) In each of the domains [A] and 
[B] the probability of finding one electron is equal, ua = 
ub = 1. (b) Interacting electrons with possible exchange, (c) 
separated electrons, and (d) several measurement domains, 

UA + nc = TlB = 1. 



\^) is a linear combination of AiJ-entangled pair states, 
W = ^ij]i^tjJ2p=ibi3\ij,P), where \ij,(3) are the Bell 
(3) states j2l| corresponding to pairs [ij] and bp are constants 

When l^f) is an eigenstate of the total spin projection 
S'jot, Eqs. l|2Jl and Q simplify further. In particular, C — 
if S^^^ = ±1, while for S^^^ — the concurrence is given 
solely with the overlap between j^*) and the particular 
AS-spin-flipped state j^") = Sj^S^^l'^), Eq. ijU, or as, 
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If probabilities for singlet and triplet are equal, the con- 
currence formula reduces to C = 2|Im^j-^-j((/3^)*(p~-| 

and if ip± = fij^''^ to C = \smS\. If the state j^-) 
corresponds to the system in continuum space, i — > 
r — (x,y,z), the only change is that summations are 
replaced by integrations of = (ri, r2; S'totl^') over 
the corresponding measurement domains, e.g., C = 
l/[A]/[s][(^+)'-(^")']d'rid\2|. 

In order to illustrate how these concurrence formu- 
lae can be applied in practice, as the first example we 
consider two interacting electrons on a one-dimensional 
lattice with iV — > oo and with the hamiltonian, Hq ~ 

To be specific, let one electron with spin "f be confined 
initially to the region A {i ^ —L) and the other electron 
in region B {i ^ L) with opposite spin, Fig. 2(a). The 
simplest initial state is two wave packets with vanishing 
momentum uncertainty Ak — > 0, the left with momen- 
tum A: > and the right with g < 0. After collision 
the electrons move apart with probability amplitude tkq 
for non-spin-flip scattering and spin-flip amplitude r^q- 
More general initial wave packets are defined with mo- 
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mentum amplitudes (pk and 4>q for spin t and |, respec- 
tively. Concurrence Eq. lO after the collision is then 
expressed as 



C 



(6) 



which simplifies to C ^ '^\tkq'rkq\ for sharp momentum 
resolution wave packets, with k = —q = k^. Note that 
C — 1 when spin-flip and non-spin-flip amplitudes coin- 
cide in accord with recent a naly sis of flying and static 
qubits entanglement E3 . 01 or of scattering of dis- 
tinguishable particles |l4| . 

Consider the prototype finite range interaction, Uij — 

\^Y^m=o^\i-3\-m- The Hubbard mod el I M = 0) can 
be solved analytically in one-dimension [25j and the am- 
plitudes are tkq = \ + rkq = (sin/c — siwq) /[{sva.k — 
sing) -f z[//(2to)]- In Fig. 2(a) concurrence is presented 
for wave packets with well defined momentum fco for 
U = to, together with a longer range interaction case, 
M = 3, for sharp momentum (full line) and for a Gaus- 
sian initial amplitude (j}k = (jj-k with Ak = tt/IO (dashed 
line). An interesting observation here is substantial re- 
duction of concurrence due to the coherent averaging in 
Eq. © . Additionally, electrons will be completely entan- 
gled at some kinetic energy comparable with the repul- 
sion, U ^ 2io(l— cos fco), where spin-flip and non-spin-flip 
amplitudes coincide. 

The concurrence formula Eq. (jSJ is derived for elec- 
tronic states when double occupancy is negligible, i.e., 
tjjjl —> 0, which in our case is strictly fulfilled only asymp- 
totically when the electrons are far apart. However, 
Eq. (jSJ can be evaluated at any time t and the result- 
ing C{t) can serve as a measure of entanglement during 
the transition from initial to final state. In Fig. 2(b) we 
present the time dependence of C(t) for some typical fcg, 
with M = 5 and U = 2tQ. Oscillation with t can be inter- 
preted as response to the finite time duration of electron- 
electron interaction and the model can be approximately 
mapped onto an effective Heisenberg model, for which 
concurrence oscillates as C{t) = | sin Jeff t|, where Jeff is 
the effective antiferromagnetic coupling between the elec- 
trons. 

Another important example is the concurrence of 
flying-static qubits in experiments in which the system 
is prepared with a static electron bound in some confin- 
ing potential (region [B]) and a flying electron injected 
in some distant region [A] [2^ . Contrary to the pre- 
vious case with translation symmetry, after the collision 
there are nonvanishing amplitudes for transmission (into 
region [C]) and reflection (back into region [A]), as shown 
in Fig. 1(d). 

Let the initial state be prepared as tp^j = {bigj ± 
gib j) 1^/2, where hi is the orbital state of the bound 
electron with spin J, centered around i ~ 0. Sim- 
ilarly, gj cx / 0fce't*''(''+-'")~'^'=*ldfc is the initial orbital 
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Figure 2: (Color online) (a) C for: (i) the Hubbard model 
(M = 0) for U = to and Afc = 0; (ii) M = 3: for Ak = (full 
line) and Afc — tt/IO (dashed), (b) C{t) for Gaussian packets 
with various fco and M — 5, U — 2to and Afc — n/20. At 
t = the separation between the packets is 2L = 10/Afc. 



state of the propagating electron with spin |, centered 
around i ~ —L and moving in the positive z-direction 
with momentum amplitude (j)k peaked at fc ~ fco, and 
with momentum uncertainty Afc — > 0. Here we con- 
sider elastic scattering with amplitudes after the colli- 
sion, (p^- = r±{biaj ±a^bj) +t±{btCj ±Ctbj), where r±(fco) 
and t±(fco) are singlet (triplet) reflection and transmis- 
sion amplitudes and aj, cj are normalised wave packets 
with mean momentum —fco and fco, respectively. 

Two basic experimental setups are possible when elec- 
trons are detected in different measurement domains, 
[AB] or [BC]. Concurrence corresponding to reflected 
qubits is then 



Cab — 



UATlB 



|r^(fco)-r^(fco)| 
\r+{koW + \r^{koW 



(7) 



where ua = {J2s ie[A] ^is), riB = 1 |2a|. Concurrence for 
transmitted qubits, Cbc, is given by an analogous ex- 
pression with A C, and consequently with r± replaced 
with t±. If the measuring apparatus captures both, re- 
flected and transmitted electrons {i € [A] U [C], j £ [B]), 
concurrence is given by Cac.b — !(''+ — +?"_) + 

(t+ — i_)(t-|- -|- i_)*| and no additional renormalisation 
is required. Eq. Q also follows directly from Eq. ^ if 
appropriately applied to scattering states [2^ . How- 
ever, for finite Afc, Cab (and correspondingly Cbc or 
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Figure 3: (Color online) Concurrence corresponding to vari- 
ous domains for infinite-?/ Anderson model with e + U — —to, 
ti ~ io/4 and Ak — > 0. Dashed dotted line represents the 
singlet transmission probability \t+\^. 



Cac,b) has to be rederived from Eq. (0), 



Cab — 



|/[r^(fc)-r^(fc)]|^fepdfc| 
/[|r+(fc)|2-t-|r_(fc)|2]|,/,,|2dfc- 



(8) 



In order to demonstrate the basic properties of Cab 
and Cbc we consider here the Anderson model, H = 
Ho +Z]s[e"Os - (^1 - to){c^_i^CQs +CosCis + h.c.)], where 
i/o is the Hubbard hamiltonian in which U — Q except for 
the impurity site, e < is the impurity energy level, and 
ti is the hopping matrix element connecting the impurity 
site 1 = with left and right leads. 

In the large-C/ regime, U,—e 3> ioi the static elec- 
tron is strongly localised, bi ^ SiQ. Electrons in the 



triplet channel are reflected, r_ 



V2' 



0, while 



singlet scattering amplitudes exhibit 'charge transfer' 
resonance: t+ = -i= + r+ = -^Tk/{ek - ujq + tTk) 

with Efc — —2t()Cosk, ll>o ~ {e + U)/{1 - 2tl/tl) and 
Tk = 2^2(4^2 _ ^2)1/2/(^2 _ 2^2) 27]. 'Transmitted' con- 
currence is due to the missing triplet amplitude, trivially, 
Cbc = 1- Reflected electrons are completely entangled 
at the singlet resonance energy but 'total' concurrence 
Cac,b = there, as shown in Fig. (0). 

The main result of this work is the closed form formulae 
of Wootters entanglement measure defined for two delo- 
calised electrons. The proposed approach enables simple 
analysis of entanglement for a variety of realistic prob- 
lems, from scattering of flying and static qubits repre- 
sented as wave packets with finite energy resolution, to 
time evolution of static qubits due to electron-electron in- 
teraction or due to externally applied fields. Further ap- 
plication to systems described with mixed states or with 
more than two electrons is possible, however, an appro- 
priate definition of entanglement valid also for systems 
with non-negliglible doubly occupancy, remains open. 
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